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ABSTRACT 

Scientific satellites may require translunar orbits aligned with the Sun- 
Earth line, with most of the period spent in either the sunward or antisun- 
ward direction. To maintain alignment, the orbit’s line of apsides must 
rotate at a rate equal to mean angular motion of the Earth about the Sun. 
To maintain this rotation of the line of apsides by use of fuel onboard the 
spacecraft is prohibitively expensive. Farquhar and Dunham (Refer- 
ence 1) proposed a method for maintaining the desired alignment by gain- 
ing momentum at the expense of the Moon during a close approach— a 
lunar swingby— as the spacecraft passes beyond lunar orbit, then return- 
ing the momentum at the second lunar swingby as the spacecraft returns 
within the lunar orbit. The cycle of double-lunar swingbys may then be 
repeated. Dunham (Reference 2) presented the orbit parameters neces- 
sary to achieve double-lunar swingby orbits which will maintain Sun-Earth 
line alignment. The details of the Keplerian approach to calculation of 
these parameters are presented. Methods for solution of the necessary 
equations for these parameters are presented. 


•This work was supported by the National Aeronautics and Space Administration (NASA) /Goddard 
Space Flight Center (GSFC), Greenbelt, Maryland, under Contract NAS 5-31500. 


1. DEFINITION OF A DOUBLE-LUNAR SWINGBY CYCLE 


Consider a Keplerian orbit about the Earth with its line of apsides pointing toward the 
Sun. After 3 months, ignoring perturbations, the line of apsides would be perpendicular 
to a line drawn from the Sun to the Earth; to maintain its Sun-pointing line of apsides, 
this line must rotate about the Earth at a rate equal to the mean angular motion of the 
Earth about the Sun. Farquhar and Dunham (Reference 1) described a method to achieve 
such an apsidal rotation rate for translunar orbits. The method requires two close lunar 
encounters, swingbys, per cycle. The first swingby occurs as the spacecraft, moving away 
from the Earth, crosses the lunar orbit with the Moon to its left. This “trailing edge” 
swingby increases the energy of the spacecraft’s orbit and hence increases the semimajor 
axis. We call the orbit with the larger semimajor axis the outer-segment loop or simply 
“outer loop,” and the original orbit the inner-segment loop or simply “inner loop.” The 
spacecraft’s outer loop orbit period is such that more than one lunar month passes follow- 
ing the first lunar swingby before the spacecraft again crosses the Moon’s orbit, this time 
moving toward the Earth. The second swingby occurs at this crossing, again with the 
Moon to the left of the spacecraft. Thus a “leading edge” swingby occurs, removing 
energy from the spacecraft’s orbit and reducing its semimajor axis to its original length— 
the length before the first lunar swingby. Next, slightly less than one lunar month passes, 
the spacecraft is now ready for another outward crossing of the Moon’s orbit, and the 
Moon and spacecraft are at the same relative position as for the first lunar swingby. This 
defines one complete “double-lunar swingby” cycle. Note that at the first swingby the 
Moon’s pull rotated the line of apsides counterclockwise. The second swingby, occurring 
as the spacecraft moved toward the Earth, also resulted in a counterclockwise rotation of 
the line of apsides. If, then, the sum of these two rotations divided by the time for one 
complete cycle equals the mean angular motion of the Sun, the spacecraft’s line of 
apsides will continue its sunward alignment. 

Figure 1 shows one complete double-lunar swingby cycle. The Moon’s positions at the 
first, second, and third (first) swingbys are shown as Si, S2, and S3. The true anomaly of 
the spacecraft at the time of the first swingby is shown as fi for the inner loop orbit and 
fo for the outer loop orbit. Thus the apsidal rotation is 2(fj - f G ) , and occurs in a time 
equal to the time from Si to S3. We can then write the first necessary condition for a 
lunar swingby as 
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Figure 1. A Double-Lunar Swingby Trajectory 
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where 

T e = Earth’s orbital period 

2t a = time spent in outer loop from Si to S 2 

t s = twice the time from perigee of inner loop to Si, plus one complete 

inner loop period 

Two more conditions relate the spacecraft’s orbital parameters to the Moon’s motion: 
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radius of Moon’s orbit 


r m = swingby distance of spacecraft from Moon 
Tm = Moon’s orbital period 

These three equations specify the geometry constraints for a complete double-lunar 
swingby cycle. 


2. REFORMULATION OF THE NECESSARY EQUATIONS 

In this Keplerian formulation, the transfer from the inner orbit segment to the outer orbit 
segment is assumed to occur instantaneously when the spacecraft crosses the lunar orbit. 
We also assume a circular lunar orbit. 
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The three equations may be expressed solely in terms of three variables; a*, Pj, and a , 
where 


a = bend angle; angle through which velocity of spacecraft is changed 
at Si 

aj = apogee distance of inner loop 

P, = perigee distance of inner loop 

The three nonlinear equations may then be solved numerically to determine these three 
unknowns. A description of the solution process is given in Section 3. We now show all 
the relations which allow Equations (1), (2), and (3) to be expressed in terms of a , a t , 
and Pj. 

For Equation (1), we require expressions for fi , fo, t a , and t s in terms of a , a i( and 

Pi- 


/«e = gravitational parameter for Earth 


ai = ^ (a; + Pj) = semimajor axis of inner loop 
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l x and y components of Moon’s velocity 
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v x = - sin f, 


v y = ei + cos f! 
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2 71 3 /2 
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For Equation (2), we express 


+ T], one inner loop period plus twice the time from 
perigee of inner loop to Si 

r m in terms of previously defined quantities. 


ju m = gravitational parameter for Moon 


Voo — {(v xs V xm ) + (Vys Vym) } 



For Equation (3), all variables have been related previously to a , a;, and Pj . 

3. SOLUTION OF THE EQUATIONS 

Writing the equations as 
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we use the Newton-Raphson method to find the required solution vector, x. 
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The method is to guess a new iterate, x N+1 , from the previous iterate, x N , by 

? n + i _ ? n _ £ (x N ) 

J(x N ) 

where J (x) is the Jacobian matrix 
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In fact, we define 

Z N+1 ■ -(x N+1 - x N ) 


and write 


J(x N ) Z N+1 = F (x N ) 

This linear system can be solved for Z N+1 once J and Fare evaluated at x N . Then x N+1 
is found from 

x N+1 = - Z N+1 + x N 

The partial derivatives are estimated numerically via central differences 


dfj = fj (xj + hjj) - f; (xj - hjj) 

9xj 2hjj 


The step size hy is chosen such that 
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dfj_ 

9xj 


fe |fil 


where f is a machine constant. For IBM double precision, e = 2.2 x 10" 16 . For IBM 
single precision, £ = 9.5 x 1(T 7 . 
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4. SAMPLE CALCULATION 


Computer output for a sample solution of a double-lunar swingby orbit is shown in the 
appendix. First, the initial guess values for the three independent variables are shown. 

a = 0.2 radians = 11.46 degrees 


aj = 700, 000 km 
Pi = 40,000 km 

Next, values for h, used to calculate partial derivatives, are displayed. Then the results of 
each iteration are shown in the form 

J (x) Z (x) = F (x) 

After seven iterations, the final values of the independent variables are shown: a = 
19.4205 degrees, a ; = 549,888 km, and Pj = 37,432 km. These results agree with those 
given on page 2-2 of Reference 2. 


c 
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APPENDIX - SAMPLE DOUBLE-LUNAR SWINGBY CALCULATION 

The following shows details for a solution of a double-lunar swingby orbit. 


PERIOD OF SATELLITE 
PERIOD OF PLANET 
PLANET GRAV. PARAMETER 
SATELLITE GRAV. PARAMETER 
SATELLITE ORBITAL RADIUS 


27.3317 DAYS , 

365.2500 DAYS , 

3 . 966E + 05 KM**3/SEC**2 , 

3.986E+05 / 81.3700) KM**3/SEC**2 

3 8<+3 9 9 . KM 


INITIAL VALUES FOR X : 

2 . OOOOOE-Ol 7. OOOOOE+05 4.00000E+04 

ERROR TOLERANCES FOR SUCCESSIVE ITERATES : 

1.00000E-06 1 . 00000E-06 1.00000E-06 


FROM HSET, H VALUES ARE: 

2 . OOOOOE-06 3 . 50000E+00 2.02500E+00 
5, 00000E-07 4.37500E-01 2.02500E+00 
6 . 75000E -06 8.75000E-01 4. OOOOOE-Ol 


***** ITERATION NUMBER 1 ***** 

THE SOLUTION TO THE FOLLOWING MATRIX EQUATION, 



3 . 0 2 2 E - 0 7 

-2.30101E-13 

2.28045E-13 | 

1 Z< X ) 

1 

1 

1 -9. 61618E-08 



-4.53911E+01 

-<*. 36530E-0 5 

| 

5. 97961E-06 | 

j 

1 Z< 2 ) 

l 

1 

1 

= 1 -2.51620E+00 



1 . 13706E +00 

1 . 29872E - 05 

2 . 72502E -05 I 

| 

1 Z< 3 ) 

1 

1 

1 

1 1.74578E+00 



-l.<f<+495E-01 2 . 0 45 96 E +05 -2.36627E+04 


HuOrV= U z^zhi E + 00 3 -^ 5E ' 01 ^ ■ 95904E + 05 6. 3 66 2 7 E+ 0^ 


***** ITERATION NUMBER 2 ***** 

THE SOLUTION TO THE FOLLOWING MATRIX EQUATION, 



5.51903E-07 

-3 . 2 93 6 6 E - 1 3 

3 . 67081E-13 | 

1 Z( 1 ) 

j 

l 

| 

1 

4 . 00435E -08 



-1 . 15960E +01 

-3 . 9 0 9 7 1 E - 05 

2 . 20171E-06 | 

1 Zi 2 ) 

1 

1 

1 

2 . 61698E+00 



7 . 5548 2 E - 0 1 

1 . 1 7588 E - 05 

2 . 18728E-05 | 

1 Zi 3 ) 

1 

1 

1 

1 

-6 . 28797E-03 



8 . 56683E -03 -6.74558E+04 3.56809E+04 


F -NORM 1 = N 8 2 89039E-01 3 ' 35878E ’ 01 s -62860E + 05 2.798X9E409 
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***** ITERATION NUMBER 3 ***** 


THE SOLUTION TO THE FOLLOWING MATRIX EQUATION, 



2 . 3 55 1 3 E - 0 7 

-5 . 67579E-13 

6 . 72816E-13 

I Z( 1 ) 

1 

i 


| -1 . 38927E-08 



-3 . 32048E+01 

-6 . 23036E-05 

1 . 39124E -05 

I Z( 2 ) 
j 

l 

l 

i 

= 

I-7.87616E-01 



4.51755E-01 

1 . 10362E-05 

2 . 66817E -05 

I Z( 3 ) 

i 

1 


1 -1 . 014 23E-01 



IS Z : -3 . 44420E-03 1.24880E+04 -8.90822E+03 


ITERATION 3, X : 3.39323E-01 5.50372E+05 3.68901E+04 

F -NORM = 5.10358E-02 


***** ITERATION NUMBER 4 ***** 


THE SOLUTION TO THE FOLLOWING MATRIX EQUATION, 


2 . 90579E-07 

-5 . 71015E-13 

5.48010E-13 1 

j 

1 Z( 1 ) 1 

j 

-4.62181E-1Q 

j 

1 -2. 78361E+01 

-5 . 74992E -05 

9.34668E-06 j 

1 Z< 2 ) 

-4 . 30255E -02 

j 5.00948E-01 

1 . 119I8E-05 

2.51765E-05 I 

1 Z ( 3 ) 1 

-8 . 01032E -03 


IS Z : 3.71385E-04 4.80824E+02 -5.39298E+02 


ITERATION 4, X : 3.38951E-01 5.49891E+05 3.74294E+04 

F-NORM = 2.36043E-04 


***** ITERATION NUMBER 5 ***** 


THE SOLUTION TO THE FOLLOWING MATRIX EQUATION, 


2 . 94573E-07 

-5. 68936E-13 

5 . 407I3E-13 

! 

i 

1 

i 

Z( 1 ) 

1 

I 

I -3 . 21530E -12 
j 

-2 . 753B2E+0I 

-5 . 71133E-05 

9 . 10926E-06 

1 

i 

f 

1 

j 

Z( 2 ) 

1 

1 

i 

= 1 -2 . 13269E-04 

j 

5 . 05424E -0 1 

I . 12043E-05 

2 . 51065E -05 

1 

1 

l 

1 

Z( 3 ) 

1 

1 -2 . 27742E -05 


IS Z : -1.49764E-08 3.35773E+00 -2.40526E+00 


ITERATION 5, X : 3.38951E-0I 5.49888E+0 5 3.74318E+04 

F-NORM = 4.30244E-09 



***** ITERATION NUMBER 6 ***** 


THE SOLUTION TO THE FOLLOHING MATRIX EQUATION, 


2 . 945 9 1 E - 0 7 

-5 . 68930E -13 

5 . 40685E -13 

1 

j 

1 

j 

Z( 1 ) 

1 

i 

I -4. 24080E-17 
1 

-2 . 75337E+01 

-5. 71117E-05 

9 . 10824E-06 

1 

1 

i 

1 

1 

1 

Z\ 2 ) 

1 

1 

j 

= I -3 .82426E-09 

j 

5 . 05440E -01 

1 . 1 2043E-05 

2.51062E-05 

1 

1 

1 

1 

Z( 3 ) 

1 

1 

I -4. 78174E-10 


IS Z : 2.49561E-11 4.83695E-05 -4.11348E-05 


ITERATION 6, X : 3.38951E-01 5.49888E+05 3.74318E+04 

F-NORM = 1 . 0 1 585 E - 1 4 


***** ITERATION NUMBER 7 ***** 


THE SOLUTION TO THE FOLLOHING MATRIX EQUATION, 


94591E -07 

-5. 68930E-13 

5.40685E-13 

1 

j 

1 

i 

Zi 1 ) 

1 

i 

1 -1 . 05879E -22 

j 

75337E +01 

-5 . 71117E-05 

9. 10824E-06 

1 

i 

1 

i 

Z( 2 ) 

1 

1 

■ 

= 1 -8 . 21565E-15 

0544 0 E - 0 1 

1 . 12043E-05 

2.51062E-05 

1 

1 

1 

1 

Zi 3 ) 

1 

1 

1 -1 . 94 2 8 9 E - 15 


IS Z : 5.44511E-17 9.81024E-11 -1.22264E-10 


ITERATION 7, X : 3.38951E-01 5.49888E+05 3.74318E+04 

F-NORM = 5. 1 6 2 54E - 1 5 

SOLUTION CONVERGED IN 7 ITERATIONS. 

X : 

3 . 38951E-01 5 . 4 9888E +05 3.743I8E+04 

Z : 

5.44511E-17 9.81024E-11 -1.22264E-10 

F : 

2.64698E-23 3.99680E-15 1.16573E-15 


ALPHA, BEND ANGLE 
APOGEE OF INNER ORBIT 
PERIGEE OF INNER ORBIT 
ECCENTRICITY OF INNER ORBIT 
APOGEE OF OUTER ORBIT 
PERIGEE OF OUTER ORBIT 
ECCENTRICITY OF OUTER ORBIT 
SHINGBY DISTANCE 


19.4205 DEGREES, 
5 . 4988764E+05 KM, 
3 . 7431801E+04 KM, 
0 .8725334 , 

8 . 9892445E+05 KM, 
1 . 0420415E+05 KM, 
0 . 7922417, 

2 . 7638929E +04 KM. 
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